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ABSTRACT 

We have calculated synthetic spectra for typical chemical element mixtures (i.e., a standard ff-enhanced distribution, and distributions 
displaying CN and ONa anticorrelations) found in the various subpopulations harboured by individual Galactic globular clusters. From 
the spectra we have determined bolometric corrections to the standard Johnson-Cousins and Stromgren filters, and finally predicted 
colours. These bolometric corrections and colour-transformations, coupled to our theoretical isochrones with the appropriate chemical 
composition, have provided us with a complete and self-consistent set of theoretical predictions for the effect of abundance variations 
on the observed cluster colour-magnitude diagrams. CNO abundance variations affect mainly wavelengths shorter than ~400 nm, 
due to the arise of molecular absorption bands in cooler atmospheres. As a consequence, colour and magnitude changes are largest 
in the blue filters, independently of using broad or intermediate bandpasses. Colour-magnitude diagrams involving uvy and UB 
filters (and their various possible colour combinations) are thus the ones best suited to infer photometrically the presence of multiple 
stellar generations in individual clusters. They are particularly sensitive to variations in the N abundance, with the largest variations 
affecting the Red Giant Branch (RGB) and lower Main Sequence (MS). BVI diagrams are expected to display multiple sequences 
only if the different populations are characterized by variations of the C+N-l-O sum and/or helium abundance, that lead to changes 
in luminosity and effective temperature, but leave the flux distribution above 400 nm practically unaffected. A variation of just the 
helium abundance, up to the level we investigate here, affects exclusively the interior structure of stars, and is largely irrelevant for 
the atmospheric structure and the resulting flux distribution in the whole wavelength range spanned by our analysis. 

Key words, stars: abundances -Hertzsprung-Russell and C-M diagrams - evolution - globular clusters: general 



1. Introduction clusters, revealed that the 'primordial' population is generally 

vo , . ... , , ^ , • , , , , N the less numerous one, comprising typically 30% of all clus- 

^ It IS now widely accepted that Galactic globular clusters (GCs) members. In contrast, 50-70% belong to the 'intermediate 
^^ host multiple stellar populations, displaying the effects of clus- j^ti^n' , defined as showing an overabundance of [Na/Fe] of 

. . ter internal chemical evolution. The majority of GCs appear ^^^^ ^^^^ 4^ -^^^ s tar-to-star error in fN a/Fe]- 

Cn , to have two approximately coeval subpopulations, or rather, determination for each cluster; see ICarretta et al.l TOl above 
^ ■ two generations of stars. The recognition of this fact came ^^^^ abundance observed in field stars with the same iron 

^ with the increasing obsei-vational evidence for what was called ^^^^^^^ ^^^^ ^^^^^^^^^ 28O8, the anoma- 

. . . star-to-star abundance variations or abundance anomalies , j.^^ ^.^^^^ ^^^^^^ ^^1^^^^ ^^.^^ ^^^^ remaining stars ai-e 

which manifest themselves as variations of nitrogen oxygen ^^^^^^^ .^^^^^^ population', and - at least in the case of 
and sodium, and in some clusters of magnesium and aluminum 28O8 (see the following discussion) - they could be repre- 

^^"^"^1"^! r*"^ ^° ^ T .'^f/ sentative of a third generation of stars, distinct from the primor- 



- 1—1 

X 

S : (NGC 6752, NG£6397, and 47 Tucanae Pasquini et al. 2005; dkfand'the [nte^medkte^onel 
- ■ ■ iBonifacioet a T 2007; Lind et al .' 2009; S hen etaIJI2010) . anti 



correlating with sodium, and hints exist that elements as heavy Such element variations in clusters have been known since a 

as sulfur might be affected (in 47 Tuc, ISbordone et al.ll2009l) . few decades (e.g., in M3 and MI3, | Cohen| | 19781) , mamly due 

These abundance variations are found in stars of all evolution- specti'oscopic analyses of cluster members, the number of 

aiy stages. An extensive discussion of these anticoiTelations, which has steadily increased, in particular due to the availabil- 

their inteipretation as two generations of cluster stai's, and the i^y of high-resolution multi-object spectrographs. Nevertheless, 

connection w i th struc tural cluster parameters can be found in *e actual fraction of cluster stars investigated is still very small. 

ICarrettaetaI.I(l20lBbl) . where also a compr ehensive list oforigi - Additional evidence that GCs host multiple populations came 

nal papers is quoted (see also the review by Gratton et al."2004'). high-qualit^otometric data for thousands of cluster 

It is remarkable that the investigation by Carretta et al. (2010b), ^t^r^- l ^ndersonl (Il998|) discovered that the broad main sequence 

which included over 1200 red giants in 19 carefully selected (MS) of the (untypical) cluster a; Cen ma^iactually be composed 

of two separate sequences. Lee et al.l (Il999h identified such sub- 

Send offprint requests to: L. Sbordone (e-mail: lsbordone@mpa- populations in the cluster red giant branch (RGB) and pointed 

garching.mpg.de) out analogies with NGC 2808 and M22 CMDs which were to be 
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confirmed almost a decade later. iBedinetalJ (12004 and follow- 
up work demonstrated that the two sequences have different 
metallicities, with the bluer one surprisingly hosting the more 
metal-rich population ([Fe/H] = -1.26 as compared to -1.57 of 
the bu lk of the cluster populatio n). The only explanation found 
so far (lNorrisl2004HPiotto et al.ll2 005) is that the bluer sequence 
has a higher helium mass fraction Y x 0.38, compared to a nor- 
mal value of 0.25. A weakly populated th ird main-se quence, red- 
der than the two others, was found by IVillanova^e t al. (2007|) 
and definit ely confirmed on the basis of an accurate WFC3 pho- 
tometry by iBellini et al.l (l2010l). ai Cen is also known for host- 
ing at least fo ur subgiant ('Lee et al. 20051 ISollima et al.ll2005t 
|yillanovaj;t al. 2007; Bellini et al. 2010) and red giant branches 
(see Johnson et al. 20091 for a recent compilation), which differ 
not only in metal content, but probably also in age. While these 
empirical findings, the cluster mass, and the undisputable pres- 
ence of a spread in iron content, definitely put ll) Cen in a special 
position among GCs, the photometric evidence for several sub- 
populations is clear, and has led to similar investigations in other 
massive clusters. 

As it turned out, multiple cluster sequences in the CMDs are 
indeed present in several other clusters. So far, they have been 
detec t ed in NGC 2808 ( three main-sequences; iD'Antona et al.l 
2005t iPiotto et alJ l200l. NGC 185 1 (double SOB and RGB 
Milone et all 120081: iHan et~an I2009 ?). M4, NGC 3201 and 



NGC 1261 (two red giant branch populations; iMarino et al.l 
l2008t iKravtsov et af] l20T0allbh . as well as in NGC 6752 and 
47 T uc (clear spre ad of the MS, also RGB spre ad in NGC 
6752 ; I Anderson et a l. 2009; Mil one et all2 010b; Kr avtsov et aP 
I20TT ). In a nu mber of cases (e.g. NGC 1261, NGC 3201, 
NGC 6752, see iKravtsov et"ani2010b.a. 2011 ). evidence exists 
that the two populations also have different radial distributions, 
with the redder RGB population being more centrally concen- 
trated. Aside from the afore mentioned NGC 185 1, [Motto (2008), 
iMarino eFal ] (12009) and M ilone et al.l (|2pi0a) announced pre- 
liminary findings of multiple subgiant branches in seven more 
clusters, among them M22, M54, 47 Tuc, NGC 6388. These 
object do not display dramatic Fe abundance spread as the one 
found in oj Cen, although M22 seems to display an abundance 
spread of about 0.2 dex, correla ting with the abundance of s- 
process elements (Da Costa et al.t2009; Marino et al..2009.) . The 
question therefore is, what the cause for the photometric split- 
ting could be? A second question concerns the fact that the split- 
ting and/or broadening of main-sequences, subgiant branches, 
and red giant branches, is not present in all photometric colors, 
a circumstantial evidence that rules out the presence of a large 
age difference between the subpopulations as the cause for the 
splitting. 

The internal chemical evolution of the clusters that led to 
the different abundances in the two stellar generations, has af- 
fected only lighter elements up to Al, and can be linked directly 
to high temperature proton-capture nucleosynthesis. In the two 
scenarios that potentially could explain the chemical properties, 
this implies an increased helium content in the material from 
which the second generation is formed. These currently dis- 
cussed options are winds from massive stars or mass loss from 
intermediate mass stars during the thermally puls ing As ymptotic 
Giant Branch phase (AGB - see, e.g.. iRenzinil l2008l and ref- 
erences therein). The increase in helium could range between 
being marginal up to the postulated value for the blue MS in 
a) Cen and the bluest MS in NGC 2808, and could be as high as 
Y = 0.38 (Piotto et al. 2007). Therefore both the spectroscopic 
and the photometric evidence for multiple cluster populations 
are intimately linked to the same mechanism. 



For these reasons it is interesting to investigate the pho- 
tometric properties of cluster isochrones as they depend on 
abundance variations. They can be used to separate subpopula- 
tions for further spectroscopic investigations, to verify known 
or suspected chemical compositions, or to assess the frac- 
tion of first and second generation stars in a cluster They 
may also help in understanding the origin of multiple MS, 
subgia nt- a nd red giant branches wi thin one cluster Salaris et al] 
(l2006h and lPietrinfemi et al.l ( l2009l) calculated stellar tracks and 
isochrones for chemical compositions showing the CN- and 
ONa-anticorrelations, as well as for enhanced helium mixtures. 
They demonstrated that in the Johnson-Cousins filters V and 
I only an extreme helium enhancement leads to a significant 
colour change as compared to a standard Pop. II mixture. They 
also predicted sUght changes in the RGB bump position, and a 
separation of first and seco r id g eneration stars along the hori- 
zontal branch. ICassisi et al.l (2008) used these isochrones to ex- 
plain the splitting of the subgiant branch (SGB) in NGC 1 85 1 as 
a cons equence of enhanced C-hN+O abundances. IVentura et al.l 
^2009*) further investigated this question, constraining the to- 
tal C-i-N-i-O abundance and the maximum helium content com- 
patible with the M S width. From the observational side, while 
Yong et al.1 (l2009t) claim a variation in this sum by a factor of 4, 



Villanova et al. [ (EOlOh dispute any evidence for this. Needless to 



say that for the AGB-pollution scenario it is of great importance 
whether the sum of the CNO-elements in the second generation 
stars is constant or has increased due to He-burning products. 



All isochrone computations mentioned so far treated the 
actual chemical composition consistently in the stellar inte- 
rior models . In particula r , app ropriate Rosseland mean opacities 
were used. I Potter et"an (1200 7) investigated in detail the influ- 
ence of individual element abundances. This limits a strict com- 
parison of isochrones to the theoretical Hertzsprung-Russell- 
diagram (HRD). However, for the precise prediction of colours, 
also stellar atmospheres and spectra with consistent chemical 
composition are needed. This is of particular importance because 
there is a clear evidence for the existence of a correlation be- 
tween splittings of evolutionary sequences in colour magnitude 
diagrams (CMDs) involving 'blue' photometric bands (such as 
the standard U Johnson band) and abundance anticorrelations, 
while these photometric peculiarities disappear when using vi- 
sual bands like V and I (see among others the case of M4, 
IMarino et al ]l2008h . 



This work presents the first determinations of bolometric 
corrections and colours for widely used blue to near infrared 
standard filters, from synthetic spectra with chemical composi- 
tions typical of both first and second generation stars found in 
GCs. Although we are postponing a detailed comparison with 
observed clusters, our results already allow a basic comparison 
with observations and specific predictions. The structure of the 
paper is as follows: The next two sections (2 and 3) present the 
methodology of our analysis, and a description of the theoretical 
model atmosphere and synthetic spectrum calculations, respec- 
tively. Application of these spectral calculations to isochrones 
representative of first and second generation stars follows in 
Sect. 4, before a summary and discussion of the results closes 
the paper. 
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2. Methodology 

We have considered a reference isochrone from the BaSTI 
databas43 (iPietrinferni et al.l|200 6'). with the following charac- 
teristics: an age of 12 Gyr, helium mass fraction Y = 0.246 
and metal mass fraction Z = 0.001, that results in an iron abun- 
dance of [Fe/H]=-1.62 for the a-enhanced metal mixture of the 
BaSTI models ([Q'/Fe]~0.4). Such an isochrone is representa- 
tive for the first generation population in a typical Galactic GC. 
A second generation of GC stars is represented by coeval (i.e. 
12 Gyr old) isochrones calculated for the same [Fe/H]=-1.62, 
Y = 0.246, but with two different choices for a metal distribu- 
tion, in which the elements C, N, O, and Na follow observed 
(anti-)correlations. We have also included an additional case for 
a mixture with CNONa anticorrelations and, additionally, an in- 
creased initial He abundance of F = 0.400 (see Fig.[l])- 

The three metal mixtures (mass fraction of metals nor- 
malized to unity) are listed in Table [1] The g-enhanced mix - 
ture employed in the BaSTI database (Pietrinfemi et al. 2006), 
which corresponds to typical first generation subpopulations in 
Galactic GCs, is labelled as 'reference'. The first mixture rep- 
resentative of second generation stars is labelled ' CNONa 1', 
and displays - compared to the reference a-enhanced mixture 
- enhancements of N and Na by 1.8 dex and 0.8 dex by mass, 
respectively, together with depletions of C and O by, respec- 
tively, 0.6 dex and 0.8 dex. This is t he same met al distribu- 
tion already used in the calculations bv lSalaris et al.l (1200 6) and 
[Pietrinferni et al. (2009). An alternative composition for second 
generation stars is labelled 'CNONa2'; it is the same as the 
CNONal mixture but for the enhancement of N that in this case 
is equal to 1.44 dex by mass. The important difference between 
CNONal and CN0Na2 'second generation' mixtures is that in 
the first case, at fixed Fe abundance, the C+N+0 mass fraction 
is enhanced by a factor of 2 compared to the reference composi- 
tion, whereas the CNONa2 mixture has the same CNO content 
(in both number and mass fractions) as the reference composi- 
tion, within 0.5%. This also implies (given that the C+N+O sum 
constitutes most of the total metal content) that for Y = 0.246 
as in the reference first generation composition, the total metal 
mass fraction Z of isochrones for the CNONal mixture has to be 
larger than the reference composition by a factor 1.84, in order 
to have [Fe/H]=-1.62. Both CNONal and CNONa2 mixtures 
are representative of extreme va lues of the CNONa anticorr ela- 
tions observed in Galactic GCs (ICarretta et al.ll2005[ l20T0bl) . In 
particular, spectroscopic measurements of the C-hNh-O sum in 
individual GCs find typically values in first and second gener a- 
tion stars that are within a factor of ~2 ( Carr etta et al. I I2OOI - 
which is the typical error bar in these estimates. 

In case of the CNONa2 mixture, the metal distribution has 
been constructed specifically to keep the same Z of the reference 
composition, when Y = 0.246. As for the CNONal metal mix- 
ture, we have also considered the case of an enhanced helium 
mass fraction Y - 0.400. Also in this case we imposed the con- 
dition of having [Fe/H]=-1.62 as in the reference composition, 
which implies a metallicity Z larger by a factor 1.46. The final 
mixtures are given in Table |2] 

Strictly speaking, our choice of keeping [Fe/H]=-1.62 for 
the CNONal mixture with enhanced helium is not in agreement 
with scenarios that could explain the chemical signatures of the 
second generation stars. In these models Fe and all other iron- 
peak elements are not affected, and indeed the absolute iron- 
abundance should be kept constant, so that an increase of Y, 

' The BaSTI database is available at 

|http : /www ■ oa- teramo . inaf . i t/BASTI | 
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Fig. 1. Theoretical isochrones from the MS to the tip of the RGB 
for the four different chemical compositions discussed in this pa- 
per The solid line corresponds to both the reference first gener- 
ation population and the second generation population with the 
CNONa2 metal mixture (both with Y - 0.246), as they are vir- 
tually identical. The black dots along this isochrone denote the 
points for which we have calculated model atmosphere and syn- 
thetic spectra. The dashed line corresponds to the second gener- 
ation population with the CNONal metal mixture (Y - 0.246), 
while the dash-dotted line is for the same population, but for an 
initial He mass fraction of F = 0.40. See the text for more de- 
tails). 



accompanied by a decrease in hydrogen abundance, therefore 
would affect [Fe/H]. So far, there is only one highly signifi- 
cant empirical evidence showing the existence of differences in 
[Fe/H] between dis tinct subpopulations w ithin a GC, e.g, the 
case of NGC 2808 jBragadiaet al .1120101) . This has been inter- 
preted as due to significant changes of the initial helium abun- 
dance between the various subpopulations, already discussed 
in the introduction. It is also important to notice that almost 
all chemical abundance studies so far were based on scaled 
solar composition model atmospheres, at most including a- 
enhancement. If a population displays a large enough departure 
from the scaled solar composition, a failure to account for the 
true chemical pattern in the model atmosphere could potentially 
skew the derived [Fe/H], possibly masking differences in [Fe/H] 
among populations. 

At any rate, even for our choice of an extreme He enhance- 
ment, by varying Y keeping Z fixed one would obtain a varia- 
tion of [F e/H] of about O.I dex. W e have verified from the a- 
enhanced ICastelh & Kuruc j (1200 4') models, that such a small 
variation of [Fe/H] affects the bolometric corrections to the pho- 
tometric filters discussed in this work at the level of at most 
~0. 01-0. 02 mag, and only at blue wavelengths. 

The HRD from the MS to the tip of the RGB of the 
isochrones considered in our analysis is shown in Fig. [1] In case 
of the CNONa2 mixture, the corresponding isochrone is prac- 
tically identical to the reference a-enhanced one, therefore we 
have used this latter isochrone to represent also a second gen- 
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Table 1. Mass and number fractions (normalized to unity) for the three metal mixtures considered. Full chemical mixtures including 
all elements from H to Ni are listed in Table |2]that appears in the online version of the paper 





Reference 




CNONal 




CN0Na2 






Number frac. 


Mass frac. 


Number frac. 


Mass frac 


Number frac. 


Mass frac. 


c 


0.108211 


0.076451 


0.013020 


0.010454 


0.027200 


0.019200 


N 


0.028462 


0.023450 


0.860012 


0.805283 


0.707661 


0.642054 





0.714945 


0.672836 


0.054256 


0.058031 


0.113300 


0.106800 


Ne 


0.071502 


0.084869 


0.034240 


0.046189 


0.071502 


0.084869 


Na 


0.000652 


0.000882 


0.001970 


0.003028 


0.0041 10 


0.005565 


Mg 


0.029125 


0.041639 


0.013947 


0.022661 


0.029125 


0.041639 


Al 


0.000900 


0.001428 


0.000431 


0.000777 


0.000900 


0.001428 


Si 


0.021591 


0.035669 


0.010339 


0.019412 


0.021591 


0.035669 


P 


0.000086 


0.000157 


0.000041 


0.000085 


0.000086 


0.000157 


s 


0.010575 


0.019942 


0.005064 


0.010853 


0.010575 


0.019942 


CI 


0.000096 


0.000201 


0.000046 


0.000109 


0.000096 


0.000201 


Ar 


0.001010 


0.002373 


0.000484 


0.001292 


0.001010 


0.002373 


K 


0.000040 


0.000092 


0.000019 


0.000050 


0.000040 


0.000092 


Ca 


0.002210 


0.005209 


0.001058 


0.002836 


0.002210 


0.005209 


Ti 


0.000137 


0.000387 


0.000066 


0.000210 


0.000137 


0.000387 


Cr 


0.000145 


0.000443 


0.000069 


0.000241 


0.000145 


0.000443 


Mn 


0.000075 


0.000242 


0.000036 


0.000132 


0.000075 


0.000242 


Fe 


0.009642 


0.031675 


0.004617 


0.017238 


0.009642 


0.031675 


Ni 


0.000595 


0.002056 


0.000285 


0.001118 


0.000595 


0.002056 



Table 2. Abundances for selected elements in the four chemical mixtures used for the atmosphere computations. The full table is 
made available in the online version. 







REFERENCE 7=0.246 


CNONal Y-- 


=0.246 




CNONal Y-- 


=0.400 




CNONa2 Y-- 


=0.246 




N 


El. 


Mass frac. 


[El]" 


Mass frac. 


[El] 


[E1]-[E1]*^^ 


Mass frac. 


[El] 


[E1]-[E1]„/ 


Mass frac. 


[El] 


[E1]-[E1]„/ 


1 


H 


7.5300e-01 


12.000 


7.5217e-01 


12.000 


0.0000 


5.9854e-01 


12.000 


0.0000 


7.5300e-01 


12.000 


0.0000 


2 


He 


2.4600e-01 


10.915 


2.4600e-01 


10.916 


0.0005 


4.0000e-01 


11.226 


0.3108 


2.4600e-01 


10.915 


0.0000 


6 


C 


7.6386e-05 


6.930 


1.9172e-05 


6.330 


-0.5999 


1.5256e-05 


6.330 


-0.5999 


1.9184e-05 


6.330 


-0.6001 


7 


N 


2.3430e-05 


6.350 


1.4768e-03 


8.150 


1.8000 


1.1752e-03 


8.150 


1.8000 


6.4150e-04 


7.787 


1.4374 


8 


O 


6.7226e-04 


7.750 


1.0642e-04 


6.950 


-0.8000 


8.4687e-05 


6.950 


-0.8000 


1.0671e-04 


6.951 


-0.7993 


11 


Na 


8.8125e-07 


4.710 


5.5531e-06 


5.510 


0.7999 


4.4189e-06 


5.510 


0.7999 


5.5602e-06 


5.510 


0.8000 


12 


Mg 


4.1603e-05 


6.360 


4.1558e-05 


6.360 


0.0000 


3.3070e-05 


6.360 


0.0000 


4.1603e-05 


6.360 


0.0000 


13 


Al 


1.4268e-06 


4.850 


1.4249e-06 


4.850 


-0.0001 


1.1339e-06 


4.850 


-0.0001 


1.4268e-06 


4.850 


0.0000 


14 


Si 


3.5638e-05 


6.230 


3.5600e-05 


6.230 


0.0000 


2.8329e-05 


6.230 


0.0000 


3.5638e-05 


6.230 


0.0000 


20 


Ca 


5.2045e-06 


5.240 


5.2009e-06 


5.240 


0.0002 


4.1387e-06 


5.240 


0.0002 


5.2045e-06 


5.240 


0.0000 


22 


Ti 


3.8667e-07 


4.034 


3.8512e-07 


4.033 


-0.0013 


3.0646e-07 


4.033 


-0.0013 


3.8667e-07 


4.034 


0.0000 


26 


Fe 


3.1648e-05 


5.880 


3.1613e-05 


5.880 


0.0000 


2.5156e-05 


5.880 


0.0000 


3.1648e-05 


5.880 


0.0000 


28 


Ni 


2.0542e-06 


4.671 


2.0503e-06 


4.670 


-0.0003 


1.6315e-06 


4.670 


-0.0003 


2.0542e-06 


4.671 


0.0000 




Z 


9.9937e-04 




1.8343e-03 






1.4597e-03 






9.9937e-04 







a [El]=logA'(El) - logA'(H) + 12 

b [El] for this mixture minus [El] for the reference mixture. 



eration population bom with the CNONa2 mixture. As we have 
teste d with specifi c calcu lations - using the same input physics 
as in ISalaris et alj ( l2006l) - as long as the C-(-N-i-0 sum is un- 
changed (in the CNONa2 mixture it is within 0.5% of the refer- 
ence metal distribution), an isochrone calculated using a metal 
mixture with typical GC CNONa anticorrelations is equal to a 
standard (T-enhanced isochrone with the same [Fe/H], which in 
case of the CNONa2 mixture corresponds to having the same 
Z. For the two CNONal mixtures with initial He abundances of 
Y = 0.248 and Y = 0.400 we used isochro nes also taken from the 
BaSTI database dPietrinferni et al.ll2009h . to represent subpopu- 
lations with extreme values of the CNONa anticorrelations and 
a large He enhancement. The CNONal, Y = 0.246 isochrone is 
essentially identical to the reference one but for the TO and sub- 
giant br anch (SGB) regions, as discussed in lSalaris et al.l (l2006l) 
and Piet rinfemi et al ] (120091) . The TO is fainter and redder, and 
the SGB is fainter compared to the reference isochrone. When 
helium is enhanced, as is well known, the MS, TO and - to a 
lesser degree - the RGB become hotter. 



Along the reference isochrone we have selected 8 key points 
(marked as black dots in Fig. [U that cover almost the full range 
of Tetf and luminosities. For each of these points we have cal- 
culated appropriate model atmospheres and synthetic spectra for 
each of the four metal mixture/He mass fraction pairs described 
before. The parameters of the model atmosphere calculations are 
reported in Table |3] The calculations denoted as 'main set' have 
been employed to produce the theoretical CMDs for first and 
second generation stars displayed in Sect. |4] while the 'test' cal- 
culations have been used to determine the eff'ect of the choice of 
micro turbulence on the results. 

The next step of our analysis has been to produce obser- 
vational CMDs in several photometric bands, starting from the 
isochrones displayed in Fig.[T] We have considered the Johnson- 
Cousins UBVI and the Stromgren uvby photometric systems, 
that cover the full wavelength range spanned by our spectra, and 
allow us to study the effect of the different metal mixtures on 
both broad- and intermediate-band filters. 

For this purpose we have first produced CMDs in se- 
lected filter combinations for all isochrones, by employing the 
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ICastelli & Kuruc d (l2004l) bolometric corre ction tables, for th e 
reference a-enhanced composition (see also lCassisi et al ]|200i- 
From our own synthetic spectra calculations, for each of the 8 
key points we have then calculated the BC differences (ABC) 
between the reference metal mixture and each of the other mix- 
tures representing second generation stars. This enables us to 
study the effect of the abundance correlations on the bolometric 
corrections at fixed T^ff and surface gravity. We have then ap- 
plied these corrections ABC - by interpolation in T^ff between 
our key points - to the isochrones with the chemical composi- 
tions of second generation stars. 

One important issue to notice is that the ABC values deter- 
mined for the chosen key points are strictly speaking not appli- 
cable to the two isochrones for the CNONal mixture in Fig. [1] In 
case of a CNONal metal mixture and 7=0.246, along the turn off 
and SGB regions, a given value of Teff corresponds to a different 
surface gravity (because of a different luminosity and slighlty 
different evolving mass) than the reference (and the CNONa2) 
isochrones. However, as we have tested with some sample cal- 
culations, these differences in gravities, that amount to 0.08 dex 
at most, do not affect appreciably the values of ABC at a fixed 

Second, all along the CNONal isochrone with Y = 0.40, the 
surface gravity differs by ~0.2-0.4 dex - at fixed Teff - from our 
key points, because of different evolving mass and luminosities. 
We have verified with some sample calculations, that changes 
in surface gravity of this order affect the corrections ABC by at 
most ~0.01 mag, and only for the u and U filters. Furthermore, 
the TO region of this He-enhanced isochrone is hotter than our 
hottest key point. We have therefore applied a linear extrapola- 
tion (by at most 150 K) to our grid of ABC values to cover the 
relevant temperature range. We have verified that we have not 
introduced any appreciable systematic error, by comparing the 
extrapolated ABC values with values obtained from appropriate 
spectra calculated for the TO point of the Y - 0.40, CNONal 
isochrone. 

Finally, we point out that a similar linear extrapolation of 
ABC (by at most 100 K) has been applied to cover the last ~0. 10- 
0.15 dex in bolometric luminosity, close to the RGB tip, a re- 
gion that is only very sparsely populated in the CMDs of typical 
Galactic GCs. To summarize, we are sure that the application of 
our ABC values to the CNONal mixture isochrones is not intro- 
ducing any significant error, which would affect our conclusions. 

3. Atmosphere models and synthetic fluxes. 

The synthetic spectra for the selected key points have been 
derived from self-consistent model atmospheres and s ynthetic 
spectra computed with the ATLAS 12 and SYNTH E (Kurucz" 
I2005al; fCastelli 2005; Sbordone 20Q1 ISbordone e t al. 2007) 
codes, respectively. ATLAS 12 employs opacity sampling to es- 
timate the line opacity in the atmosphere, allowing the com- 
putation of model atmospheres for arbitrary chemical compo- 
sitions. For each considered chemical mixture and each set 
of atmospheric parameters (Teff, iogg, Vturb) an ATLAS 12 
plane-parallel, LTE model atmosphere was computed, and a 
SYNTHE spectral synthesis was performed between 300 nm 
and 1000 nm. The spectral synthesis acco unts for the full s et of 
atomic and molecular lines included in the'Sbo rdond(l200l and 
[sbordone et al. (2007) Linux ATLAS porj^ as well as the pre- 



dicted lines no rmally used for pre-tabulated opacities calculation 
(lKuruczl l2005b). for the purpose of this spectral synthesis was 
to derive intermediate-to-broad band colours. All the computed 
models and synthetic spectra will be made available online to 
the community 

3. 1 . Producing the full chemical mixtures 

The first step of our calculations was to provide the full chemical 
mixtures for the atmospheres. The three chemical compositions 
of Table [T] were derived from the stellar interior calculations, 
and included only those elements relevant for that purpose in 
the form of mass fractions of metals normalized to unity. We 
will henceforth use the term 'interior' when speaking about this 
restricted set of elements. 

While it is possible to consider only a subset of elements for 
the stellar interior modeling without adversely influencing the 
results, abundances for all remaining metals should be included 
in synthetic spectrum calculations. The fractions of elements in 
the interior mixture should then be converted into absolute abun- 
dances, including appropriate values for H and He, and the miss- 
ing elements added. This process is essentially arbitrary, since 
there is no single way to fill in the missing metal abundances. 
We proceeded as follows: 

- The interior metal mixture, as listed in Table [T] was con- 
verted to a full metal mixture by filling in all the elements 
except H and He. To do so, we took as reference the chem- 
ical mixture of an g-enha nced, JFe/H1=-l.5 model from the 
ICastelli & Kuruczl (|2004|) gricQ The total mass fraction of 
the elements considered for the interiors amounts to 0.986 
in this mixture. Thus, their abundances in the interior mix- 
tures of Table [1] were scaled by this factor, and the reman- 
ing 1 .4% distributed over the missing elements according to 
their abundance fractions in the Castelli mixture. This has 
produced the three full metal mixtures of Table |2] 

- The helium mass fraction Y was set to the two possible val- 
ues, Y = 0.246 or Y = 0.400. To produce a full set of abun- 
dances we then only needed to set a ratio between the H mass 
fraction X and the cumulative mass fraction of metals, Z. 
This was done by enforcing [Fe/H]=-1.62. We came to such 
choice since observational evidence is that stars from differ- 
ent populations in GCs usually share the same [Fe/H] (see 
also the discussion in Sect. 2). The final full mixtures used 
for the atmospheres are listed in Table |2] 

3.2. The model atmospheres 

ATLAS 12 has been employed to produce sets of model atmo- 
spheres from lower MS to bright RGB stars - the parameters are 
listed in Table |3]- for each of the four chemical compositions 
given in Table |2] 

Figure |2] compares the temperature stratification in models 
employing the reference and CNONal (Y - 0.246) mixtures. 
Shown are temperature (left panel) and the corresponding rel- 
ative differences (right panel) against Rosseland optical depth 
for three models corresponding to a cool MS, a TO and an upper 
RGB star Significant differences arise only in the outermost lay- 
ers and only in cool models. As we will show below, these differ- 
ences arise from the generally lower molecular absorption in the 



^ We did not include TiO lines in the calculations, after verifying that 
no significant TiO bands formed in some representative cases, given 
that they would have made the calculations remarkably heavier. 



^ For this preprint, we have not yet set up a download system which 
will be made available for the final publication. In the meantime, anyone 
interested might contact us directly. 

Available online at http : //wwwuser . oat ■ ts ■ astro . it/castelli/gri( 



6 



Sbordone et al.: Multiple stellar populations in globular clusters 




Fig. 2. The left panel displays the temperature structures for a MS, TO, and a RGB atmosphere as used in the present work. The 
black continuous line represents the model computed with the reference chemical mixture, the red dashed line the one assuming the 
CNONal mixture. In both cases, Y - 0.246. The right panel shows the percentage temperature difference between the two different 
chemical compositions. 



Table 3. Parameters of the model atmospheres. Calculations 
with the 'main set' of parameters have been employed to pro- 
duce theoretical CMDs for first and second generation stars. 
Calculations with the 'test' parameter sets have been used to es- 
timate the effect of microturbulence. 







Vturb 




Mixtures 




K 


cm/s- 


km/s 








main set 


4100 


0.50 


2.0 




all 




4476 


1.20 


2.0 




all 




4892 


2.06 


2.0 




all 




5312 


3.21 


2.0 




all 




5854 


3.78 


2.0 




all 




6490 


4.22 


2.0 




all 




6131 


4.50 


2.0 




all 




4621 


4.77 


2.0 




all 




Vtuib test 


4476 


1.20 


0.5 


ref. Y= 


=0.246, CNONa2 Y= 


=0.4 


4476 


1.20 


1.0 


ref. Y= 


=0.246, CNONa2 Y= 


=0.4 


4476 


1.20 


1.5 


ref. Y= 
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=0.4 


6490 


4.22 


0.5 


ref. Y= 


=0.246, CN0Na2 Y= 


=0.4 


6490 


4.22 


1.0 


ref. Y= 


=0.246, CNONa2 Y= 


=0.4 


6490 


4.22 


1.5 


ref Y= 


=0.246, CNONa2 Y= 


=0.4 


4621 


4.77 


0.5 


ref Y= 


=0.246, CN0Na2 Y= 


=0.4 


4621 


4.77 


1.0 


ref Y= 


=0.246, CNONa2 Y= 


=0.4 


4621 


4.77 


1.5 


ref Y= 


=0.246, CN0Na2 Y= 


=0.4 



reference mixture compared to the other models, which affects 
mostly the outer layers of the cooler atmospheres. The strongest 
differences arise in the cool giants and are still appreciable in the 
coolest dwarf models, but become negligible in the warmer stars 
close to the TO or along the SGB. In the Teg - 4476 K giant, 
differences range from about 50 K at Tross^lO"'* to about 250 K 
at Tio,ss=10"^, which is the highest value encountered amongst 
the models explored in this work. 

Figure [3] displays the temperature stratification of the RGB 
models shown in Fig.[2]for all our four mixtures, as well as their 




-6 -4 -2 



Fig. 3. The upper panel compares the temperature stratification 
for the four Teff - 4476 K, log g= 1 .2 giant models. The solid line 
denotes the reference mixture, the dashed line the CNONa mix- 
ture with Y = 0.246, the dot-dashed Une the CNONal mixture 
with Y - 0.4, and the dash-dot-dot line the CN0Na2 mixture. 
In the lower panel the percentage temperature difference for the 
last three mixtures with respect to the reference one is plotted; 
colour and line type are coded as in the upper panel. 

deviations from the case of the reference mixture. The results 
are representative of the general result for all models. Even large 
variations in the He abundance have a quite minor effect on the 
atmosphere structure (as is evident from the extremely similar 
structures in the CNONal Y - 0.246 and Y - 0.4 cases). It 
might be worth reminding that this does not mean that the he- 
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Fig. 4. The flux distribution for the RGB model with Tejf - 4476 K and log g-1.2 for the reference mixture (black) and the CNONa 1 
mixture (red), both with the standard helium content of Y - 0.246. We have superimposed the transmission curves for the Johnson- 
Cousins U, B, V, and I filters (thin, black lines; left to right) and for the Strogren uvby filtes (grey-shaded regions). The synthetic 
spectra are broadened by convolution with a gaussian of FWHM=1700 km/s for readability purposes, and in order to roughly 
match the resolution of the Castelli & Kurucz (2004) ATLAS9-based fluxes. A number of molecular bands which vary significantly 
between the two mixtures are labeled by the name of the corresponding molecule. 



lium content will not affect the colours of second generation 
stars, since varying the He content has an important effect on 
the stellar interiors (Fig. [T). It merely indicates that the helium 
content does not affect these cool atmospheres, and thus the re- 
lationship between stellar parameters and colours. Passing from 
the CNONa 1 to the CNONa2 mixture has a slightly more no- 
ticeable effect, but we want to stress that the difference never 
exceeds 25 K. 

As far as the model atmospheres are concerned, one can thus 
conclude that the explored chemical mixtures do not have a ma- 
jor effect on the atmospheric structure, except in the outermost 
layers, where - at any rate - LTE models such as the ones em- 
ployed here are less physically sound. 

3.3. The synthetic spectra 

We display in Figs.|4]|5] and|6]the synthetic fluxes computed for 
the same three sets of parameters as in Fig. |2] as representative of 
the effects of changing chemical composition. More specifically, 
we display the case of the reference (black) and CNONa 1 (red) 
mixtures, both with Y = 0.246. Synthetic spectra were com- 
puted with very fine sampling (AA/A = 300000) but are shown 
here after convolution with a FWHM=1700 km/s Gaussian, for 
re adibility purposes, and in o rder to roughly match the resolution 
of ICastefli & Kuruc3 (l2004l) fluxes. In each plot we label some 
prominent absorption features (mostly molecular bands), whose 
strength changes significantly among chemical mixtures. 

Figure |4]compares the synthetic spectra of a Tetj - 4476 K, 
log §=1.2 giant. The CNONal mixture shows much stronger NH 
and CN absorption bands in the U, B, and I filters. This is es- 
sentially due to the much higher N abundance of the CNONal 
mixture, despite the fact that the C abundance is lower This in- 
dicates that the N abundance acts as bottleneck in forming CN 



molecules. Conversely, the G-band, a CH feature falling into the 
B filter, appears stronger in the reference mixture, since the C 
abundance is higher The increased opacity in the blue part of 
the spectrum for the CNONal mixture leads to the increase of 
the continuum flux redwards of about 450 nm, which explains 
the overall higher flux observed in the CNONal spectrum be- 
tween 450 and 690 nm, and further in the red in the intervals 
along the strong CN absorption bands. 

In the spectrum of the TO atmosphere (Teff - 6490 K, 
log §=4.22, Fig. |5]i the much higher temperature prevents the 
formation of CN and CH alike, thus removing the most promi- 
nent causes for increased blue opacity in the CNONal case, as 
well as most of the influence of the G-band. Only the very strong 
NH band around 340 nm is still visible. The much more similar 
spectra for the two mixtures in the TO atmosphere are reflected 
in the extremely similar atmospheric structures displayed by the 
two TO models in Fig. |2] 

The cool MS star (Teg = 4621 K, log§=4.77. Fig. |6l) dis- 
plays features similar to the RGB one, but shows also a much 
stronger OH absorption at the blue edge of the U filter range in 
the reference case, due to the higher O abundance: the same band 
is visible, but much less prominent, also in the MS and RGB 
stars. The Na D doublet becomes strongly wing-dominated in 
this spectrum and thus appears much stronger in the CNONal 
spectrum, where the Na abundance has been increased to repro- 
duce the observed Na-O anticorrelation. Red CN bands have an 
almost negligible effect in this star. The strong absorption feature 
around 415 nm is the MgH A-X band. 

In Fig. |7] we finally show how the blue-visible part of the 
spectrum changes among all the four mixtures considered, in 
the same RGB stellar model as described above. Two effects are 
immediately apparent: first, the CNONa2 mixture - sharing the 
same Z as the reference one, as well as a lower N enhancement - 
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Fig. 5. As in Fig.|4]but now for the TO model with Tes = 6490 K and logg=4.22. 
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Fig. 6. As in Fig.|4]but now for the TO model with T^ff = 4621 K and logg=4.47. Here, the Na D doublet is also labeled 



shows less prominent NH and CN bands than the CNONal mix- 
ture, as well as a continuum which resembles more the one of the 
reference mixture. Second, the change in the He abundance be- 
tween the two CNONal mixture has only a tiny effect on the flux 
distribution. It has here to be considered that our choice of con- 
serving [Fe/H] among mixtures leads to a lower Fe abundance 
(and Z, for that matter) in the CNONal mixture with Y - 0.400 
with respect to the one with Y - 0.246, thus mitigating the main 
effect of increasing He, i.e. an increase in molecular weight. 

Our synth etic fluxes have been co mpared with MARCS syn- 
thetic fluxes ("Gustafss on et al.l 12008). obtaining equivalent re- 
sults for the reference composition. MARCS "CN-cycled" fluxes 
are produced with CNO abundances that depart from the stan- 
dard mixture less than the ones we considered here, which, as 



expected, leads to a behavior which is intermediate between our 
reference and CNONal/2 cases. 

4. Bolometric corrections and 
colour-magnitude-diagrams 

As discussed in Section |2] we have calculated the BC dif- 
ferences ABC for each of the key points, in the UBVI, and 
uvby photometric filters. We followed the method pres ented in 
Girardi et al. (200 2|) using passband definitions from iBesselll 
( 1990) and Stromgren (1956). The results for the different chem- 
ical compositions are displayed in Figs. [SlfTOl Figure |8] shows 
ABC for the CNONa2 mixture as a function of logg. As ex- 
pected from the comparisons of the previous section, the U and 
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Fig. 7. The range between 300 and 600 nm is displayed for all 
mixtures considered, in the case of the Tgff = 4476 K, log ^=1.2 
giant atmosphere. Black line: reference mixture; red: CNONal 
(Y = 0.246); blue: CNONal (F = 0.400); green: CNONa2 (F = 
0.246). 



u filters are the most affected by the different metal mixtures. 
The BC for the reference mixture is systematically higher by an 
amount that depends on the evolutionary phase. In both U and u 
bands ABC decreases moving from the low main MS towards the 
TO (where it is practically negligible) then increases again along 
the RGB, reaching a local maximum around Teff - 4476 K and 
log g= 1.2, before decreasing slowly towards the tip of the RGB. 
Here, the largest values of ABC are of the order of 0.2 mag in u 
and 0. 1 mag in U. The effect on the other filters is much smaller: 
the largest values of ABC are attained along the bright RGB and 
are of the order of 0.02-0.04 mag. Figure |9] displays the same 
comparison for the case of the CNONal mixture, the mixture 
with enhanced C+N+O. The behaviour of ABC is qualitatively 
identical to the previous case, but the maximum values of ABC 
are increased. They are equal to -0.30-0. 35 mag in u, ~0.2 mag 
in U, ~0. 1 mag in v and ~0.04 mag in the remaining filters. The 
final comparison among BCs is displayed in Fig.[TOl that shows 
the differences between the CNONal mixture with Y = 0.246, 
and the same metal mixture with a much higher Y = 0.40. The 
values of ABC are reduced compared to the case of varying the 
metal mixture. ABC is essentially zero in all filters from the low 
MS to the lower RGB, then becomes more and more negative 
while climbing the RGB. A higher helium abundance appears to 
increase the value of BC at low Tetf and low gravities. The largest 
differences are obtained for the U and u bands, but are of the or- 
der of only 0.05 mag at most. These variations are small with re- 
spect to the effect of the metal mixture in U and u, but much more 
comparable for the other filters. Overall, th e results of Fig .[TOlar e 
broadly consistent with the conclusions bv lGirardi et al.l (l2007h . 
who studied the effect of enhanced He (for a scaled solar metal 
mixture) on the BCs for the UBVRIJHK filters. 

Finally, we have also investigated whether the values of ABC 
are affected by the choice of microturbulence, using the addi- 
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Fig. 8. Upper panel: Difference between the UB VI bolometric 
corrections ABC for the reference mixture and the CNONa2 
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panel, but for the uvby system. 
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Fig. 9. As in Fig.|8]but for the CNONal mixture. 



tional 'test' calculations listed in Table |3] Microturbulence af- 
fects line absorption by reducing line saturation: the net effect 
is that saturated lines absorb more flux when microturbulence is 
higher, while unsaturated lines do not change. Since typical mi- 
croturbulence values vary with the stellar evolutionary phases, 
one might need to take it into account in estimating ABC. In 
practice, we find that such effect is negligible: changing the value 
of Viurb from 2 km/s to 0.5 km/s changes the values of ABC by 
at most 0.01 mag. 
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Following the method outlined in Sect. |2] we have then per- 
formed the first fully consistent theoretical study of the effect 
of the CNONa abundance anticoiTelations on the MS-TO-RGB 
CMDs of Galactic GCs. In particular, we are now in the position 
to account for the CNONa pattern both in interior models and 
bolometric corrections. 

Figures [TT] to [15] display several CMDs in Johnson-Cousins 
as well as in Stromgren filters, for the isochrones of Fig. [1] Given 
that the CNONa pattern we considered for the second generation 
stars is characterised by extreme values for the anticorrelation 
observed in Galactic GCs, the range of colours spanned by our 
isochrones should give a rough idea - considering also that the 
extension of th e abundance anticorr elations varies from cluster 
to cluster (e.g. ICarretta et al.ll2010bl) - of the maximum colour 
spread to be expected in the CMD of a generic Galactic GC. 

The commonly used BVI diagrams are shown in Fig.[TT] The 
behaviour of the isochrones in these CMDs mirrors closely the 
HRD of Fig.[T] the reason being that in these filters the BCs are 
hardly affected by the change in the metal mixture and initial He 
content. As long as the sum of C-hN-i-O is constant, second gen- 
eration stars are expected to overlap with first generation objects. 
In case of enhanced C+N+O, only the TO and SGB regions are 
affected, whereas MS and RGB remain unchanged. An increase 
of Y shifts the MS (and to a lesser degree the RGB) towards bluer 
colours because of hotter Teff in the evolutionary models. 

The situation is - not unexpectedly - different when con- 
sidering the UBV CMDs displayed in Fig. [12] given that U is 
the filter most affected by the change of the metal mixture due 
to the emergence of strong molecular absorption in the atmo- 
sphere. In both diagrams the four isochrones are well separated 
along the various branches. Overall, our isochrone representative 
of second generation stars with enhanced C-hN-hO is the red- 
dest (dash-dotted line). Notice how second generation stars bom 
from a metal mixture with constant C-i-N+0 (long-dashed line) 
follow a distinct sequence from the first generation population 
(solid line), in spite of the fact that the HRDs of the underlying 
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in Fig. [T] Different line-styles correspond to different element 
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isochrones are identical. This is entirely due to the effect of the 
metal mixture on the BCs. 

The largest differences among the various isochrones appear 
along the RGB, where the mixture with anticorrelations causes 
redder (U-B) and (U-V) colours, at fixed M\j. For example, at 
Mu=2.0 the RGBs representative of second generation stars are 
redder by up to ~0.2 mag in (U-B) and ~0.3 mag in (U-V), 
depending on the metal mixture considered. We will compare 
this result with observational evidence in the discussion section. 

It is also important to notice that an increase of Y up to 0.40 
- a very extreme He-enhancement - shifts the isochrones of sec- 
ond generation objects, which tend to be redder than the ref- 
erence isochrone due to the CNONa-variations, bluewards and 
thus closer again to the reference isochrone. This is an effect 
of the stellar interior models only, and may potentially even pro- 
duce a bluer MS - depending on whether the anticorrelations are 
accompanied by a C+N+O enhancement or whether the C+N+O 
sum stays unchanged - especially in the U-(U-V) diagrams. The 
effect can best be seen by comparing the dot-dashed and dotted 
line in Fig. [12] 

As for the Stromgren filters, the M,.-(m - y) and My-(v - y) 
CMDs are displayed in Fig. [13] The CMD with (v - y) colours 
show a behaviour very similar to the case of the V-(B-V) plane, 
except for the RGB. When the initial Y is kept constant, the 
bright RGB for the C+N+O enhanced mixture (dash-dotted line) 
is separated from the reference isochrone (solid) - at My = -2.0 
the difference in (v - y) is ~0.2 mag - whereas the mixture 
with anticorrelations at constant C+N+O (long-dashed) does not 
modify the colour of the RGB. This difference in the behaviour 
of the RGB (v - y) colour for the two different mixtures with 
CNONa anticorrelations is due essentially to the different BCs, 
rather than to the underlying isochrones. 

Finally, the My-(u - y) CMD is similar qualitatively to the 
case of the U-(U-B) diagram, with all different sequences well 
separated in the CMD when the initial He is kept constant. 
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(U-B) (U-V) 

Fig. 12. As in Fig. [lUbut for the Mu-(U-B) and Mu-(U-V) 
CMDs. 
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Fig. 13. As in Fig.fTTlbut for the My-{v-y) and My-{u-y) CMDs. 



Figure fT4l displays the m\-{u - y) diagram, where the mea- 
sured m\ - (v -h) - {h - y)v!\ RGB stars is often used as an es- 
timator of the cluster metallicity (see, e.g.. 'Calamida et al.l l2007i 
and references therein). Our comparison shows that the presence 
of CNONa anticorrelations produces a spread in the location of 
GC stars in this diagram. Focusing our attention on RGB stars 
only (the right part of the curves), the four different isochrones 
in Fig.[T4]display quite different RGBs, which intersect at {u-y) 
around 2.6-2.8 mag. It is very interesting to notice that even in 
case of a constant C-i-Nh-O sum, [Fe/H], Y and total metallicity 
Z (the CNONa2 case; long-dashed line), the RGB representative 




base RGB 



2.0 3.0 4.0 

(u-y) 

Fig. 14. As in Fig.Obut for the oti-(m - y) CMD. 

of second generation stars is different from the reference RGB, 
a clear effect of the difference in BCs. 

Finally, we consider the behaviour of the index Cy-c\ -(b-y), 
where ci - (m - v) - (v - b) is found em pirically to be sensitive to 
the N abundance (see, e.g. jYong et al.l 2Q08). The c,, index repre- 
sents well ci, but removes much of the temperature sensitivity of 
this index. As a result, empirical V-c,. diagrams of Galactic GCs 
display an a lmost vertical R GB at luminosities lower than the 
RGB bump (" Yong et al]|2008 ). Figure[T5]shows how this empir- 
ical behaviour is reproduced by the theoretical isochrones. The 
theoretical RGBs of the isochrones with our four selected chem- 
ical compositions are approximately vertical below the bump lu- 
minosity, and are very well separated in c, colour All isochrones 
representative of second generation stars have redder Cy colours, 
and the isochrone with the largest N abundance (at fixed Y - the 
dot-dashed isochrone for the CNO Na 1 mixture) is the reddest, in 
agree ment with empirical results (iLind et al.ll2010l lYong et al] 
l2008l) . Notice how in this diagram an increase of Y tends to 
move the RGBs of isochrones with CNONa anticorrelations fur- 
ther away from the reference isochrone (dotted line), despite the 
fact that the CNONa 1 Y=0.4 mixture has a lower N abundance 
than the CNONal Y=0.246 mixture. 



5. Discussion 

Using ATLAS 12 model atmospheres and SYNTHE spectrosyn- 
theses we calculated synthetic spectra for typical mixtures found 
in the various subpopulations of globular clusters. From the the- 
oretical spectra we have determined bolometric corrections for 
standard Johnson-Cousins and Stromgren filters, and finally pre- 
dicted colours. 

As a result of this work, we now have in hand complete and 
self-consistent theoretical predictions of the effect of abundance 
variations on observed CMDs, as they are found or discussed 
for the different globular cluster stellar generations. Globally, el- 
ement abundance variations affect mainly the part of the spec- 
tra short of about 400 nm due to changes in molecular bands. 
Therefore, colour changes are largest in the blue filters and their 
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Fig. 15. As in Fig.Obut for the Mv-(cv) CMD. 

detection is easiest in the blue, independent of using broadband 
or narrow filters. Enhanced helium abundance is affecting only 
the interior structure of stars, leading to changes in luminosity 
and effective temperature, but is irrelevant for the atmospheric 
structure at fixed log g and T^ff, even for Y as high as 0.400. 

In this study, we did not investigate what the effect of CNO 
and He abundance would be in cluster with different [Fe/H]. 
While we don't expect the general behavior to change, the quan- 
titative result is bound to be sensitive to the overall metal content 
of the cluster 

Using the combination of the previously computed 
isochrones and the new bolometric corrections for the different 
GC subpopulations, we summarize here the effects on the var- 
ious possible colours in both Johnson-Cousins and Stromgren 
filter systems: 

1. BVI-diagrams: 

- a splitting of sequences along the MS up to the TO, and 
to a lesser degree of the RGB can only be achieved - 
within the element variations discussed in this paper - 
by varying the helium content Y. The CNONa anticorre- 
lations influence neither the stellar models nor the spec- 
trum sufficiently when the C-nN-i-O-abundance is un- 
changed (mixture CNONa2 - Fig.fTTTi: 

- on the other hand, a variation of the C+N+O-abundance 
(mixture CNONal) leads to a split of the SGB; this is 
entirely an effect of the stellar models. 

2. UBV- and uy-diagrams: anticorrelations in CNONa abun- 
dances as well as ^-differences may lead to multiple se- 
quences from the MS to the RGB, where the effect tends 
to be larger, and may reach 0.2-0.3 mag. This multiplicity is 
independent of the sum of C+N+O (Fig.[T2j left panel). The 
individual element variations are decisive. Helium enhance- 
ment, however, works in the opposite direction than CNONa 
anticorrelations (Fig.[T3j left panel). 

3. vy -diagrams: 

^ In real globular clusters the splitting may in fact appear as a spread, 
depending on the nature of the element abundance variations 



- as in the case of the fiVZ-colours, a splitting of the MS 
up to the TO can be achieved only by a variation in Y; 

- similarly, after the TO, a split of the SGB is the result of 
a change in C-i-N+O; 

- additionally, a split along the RGB may result both from 
helium and from C+N-hO variations; this is different 
from the Z?y/-case (Fig.[T3j right panel). 

4. mi uy-diagrams: 

- CNONa anticorrelations lead to splits along the MS; 

- along the SG and RGB the same anticorrelations, but also 
helium variations lead to colour differences; 

- the sign of the colour change is different for the lower 
and upper part of the RGB (Fig. HM . 

5. Cy V-diagrams: here, all parts of a CMD show the influence of 
both element anticorrelations and of helium variations, and a 
strong separation can be seen (Fig. [15). 

Our calculations also enable us to br iefly address the issue 
raised by the photometric observations bv lLee et al.l(l2009l) . who 
find broad distributions of the hk index - defined as hk-{Ca-b)- 
(b - y)- in RGB sequences of several Galactic GCs; they deter- 
mined a broadening of the order of ~0.2 mag in hk at fixed V 
magnitude and attribute this to a spread in calcium abundance, 
given that the hk index is sensitive to the ionized calcium H 
and K lines. On the other hand, direct spectroscopic measure- 
ments of calcium abundance by Carre tta et al. (2010a) do not 
find any significant spread in individual clusters. We have tested 
with our spectra whether the CNONa abundance anticorrela- 
tions affect the hk index at our chosen RGB key points. For this 
purpose we have calculated differences in the hk index among 
our adopted chemical compositions employing the appropriate 
profile fr om the Ca filter taken fr om the Asiago photometric 
database (iFiorucci & Munarill200l . We find an effect of at most 
only ~0.04 mag in the hk colour at fixed My, second generation 
stars being systematically bluer Therefore, CNONa anticorrela- 
tions cannot fully explain /i^-variations of 0.2 mag. 

In comparison with observations we find that our My - Cy- 
diagram closely follows empirical results by 'Yong e t aP (l2008h 
and Lind et al. (2010), which adds support to our theoretical 
models. We thereby also confirm that this index is well suited 
to investigate photometrically variations of nitrogen abundance 
along the lower RGB. 

Also the U-(U-B) CMD of Fig.[T2]interestingly agrees with 
the results by Marino et al. (2008) about M4 RGB stars. These 
authors have found that in the U-(U-B) CMD, objects with high 
Na - represented by our isochrones for the two compositions 
with CNONa anticorrelations - are distributed to the red of stars 
with low Na - represented by our isochrone with the reference 
chemical composition. 

A further interesting cas e is given by NGC6752. In Sect. 1 
we already mentioned that iMilone et al.l (1201 Obi) report about a 
MS splitting in visual HST filters, and an RGB spl it in a UB- 
CMD. The SGB appears to be single in the visual. lYong et al] 
(2008) found a variation in nitrogen abun dance of up to 1 .95 dex , 
correlating with the ci index. Recently, iKravtsov et alj (1201 lb 
claimed a splitting of the SGB in the U vs. B-I diagram, and a 
broadening of the RGB in the U vs. U-B plane. The single SGB 
in B-I is indicative of a constant C-i-N+O, the MS splitting of 
He-enhancement in the second generation stars, and the broad- 
ening of the RGB in blue filters agrees with both. In conclusions, 
if these findings are confirmed, they would indicate a subpopu- 
lation of CN0Na2 type, but with enhanced helium. 

Finally, we comment on the case of NGC 1851 that has been 
observed in different bands. The y-(v-y), y-(b-y), and mi-(u-y) 
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CMDs by ICalamida et alJ (l2007h show a spread, if not a bi- 
modality, along the RGB. The HST otf606W-(wf6()6w - wf814w) 
diagram by Milone et al. (2008), approximately equivalent to 
our Mv vs. V-I CMD, displays a d ouble SGB but no other spUt- 
ting. In the UBVI observations o f iHan et al.1 (|2009|) . finally, no 
RGB split is observed in the My vs. V-I CMD, whil e it is quite 
evident in the Mu vs. (U-I) one (see ' Han et"ani2009 . Fig. 1). It 
has to be noted that lHan et a l. (2009) could not detect any SGB 
split in the My vs. V-I CMD, probably due to the lower quality of 
th eir ground-ba sed photometry compared to the HST-based one 
in iMilone et al.l (|2008). On the basis of our results, only a sec- 
ond generation population with enhanced C+N-i-O and 'normal' 
He can satisfy all these constraints. In fact, a second generation 
with constant Ch-N-hO will not produce a spread along the RGB 
in the vy CMD, and a He enhancement would reveal itself along 
the MS of the »iFfin6w-(»i Ffinfiw - w^fsi^w) CM D. Also the syn- 
thetic HB simulations bv iSalaris et al. (2008) seem to exclude a 
He variation among the cluster subpopulations. 

In summary, we are now able to investigate in detail clusters 
with known or suspected abundance variations and several sub- 
populations, using information from their CMDs, on the basis of 
fully self-consistent interior and atmosphere models. As a first 
application, in a following paper we will analyze in detail the 
cluster NGC 1851 with more specific calculations. 
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a [El]-log7V(El) - log/V(H) + 12 

b [Ell for this mixture minus [El] for the reference mixture. 



